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Abstract 

We use QCD sum rules with nonlocal condensates to re-calculate more accurately the moments 
and their confidence intervals of the twist- 2 pion distribution amplitude including radiative cor- 
rections. We are thus able to construct an admissible set of pion distribution amplitudes which 
define a reliability region in the ai, 04 plane of the Gegenbauer polynomial expansion coefficients. 
We emphasize that models like that of Chernyak and Zhitnitsky, as well as the asymptotic so- 
lution, are excluded from this set. We show that the determined 02, 04 region strongly overlaps 
with that extracted from the CLEO data by Schmedding and Yakovlev and that this region is 
also not far from the results of the first direct measurement of the pion valence quark momentum 
distribution by the Fermilab E791 collaboration. Comparisons with recent lattice calculations and 
instanton-based models are briefly discussed. 
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I. INTRODUCTION 



The pion distribution amplitude (DA) of twist-2, ^(x, /x 2 ), i.e., the integral over trans- 
verse momenta of the renormalized valence-quark wave function on the light cone, is a gauge- 
and process-independent characteristic of the pion and, due to factorization theorems 0,0, 
it enters as the central input various QCD calculations of hard exclusive processes. A re- 
liable derivation of the pion DA from first principles in QCD is therefore an outstanding 
problem of paramount importance, given that it contains all of the bound-state dynamics 
and specifies in a universal way the longitudinal momentum xP distribution of the valence 
quarks in the pion with momentum P (see, e.g., for a review), 



(0 | d(z)y»y B u(0) | tt(P)> = if v P» 1 1 dxe lx(zP ^ (p„(x, /i 2 ) . (1.1) 







Recently, the CLEO collaboration [4] has measured the 7*7 — > 7r° form factor with high 
precision. These data sets have been processed by Schmedding and Yakovlev (S&Y) [5J] 
using light-cone QCD sum rules and including perturbative QCD contributions in the NLO 
approximation (NLA), to obtain useful constraints on the shape of the pion DA in terms 
of confidence regions for the Gegenbauer coefficients a 2 and 0,4. On the other hand, the 
first direct measurement of the transverse momentum distribution in the pion via diffractive 
dissociation into di-jets by the Fermilab E791 Collaboration supplemented by new lattice 
results 0] for the same quantity, also provides the possibility to deduce the shape of the 
pion DA (modulo inherent method uncertainties). 

These important findings inevitably raise the question of whether pion DAs, derived from 
first principles of QCD, can confront them successfully. The present analysis is a targeted 
investigation of these issues, the focus being on pion DAs, reconstructed from QCD sum 
rules with nonlocal condensates (NLC-SR). 

This approach, developed in fl, 0, [l(], EH by A. Radyushkin and two of us (A.B. and 
S.M.) provides a reliable method to construct hadron DAs that inherently accounts for the 
fact that quarks and gluons can flow through the QCD vacuum with non-zero momentum 
k q . This means, in particular, that the average virtuality of vacuum quarks (A; 2 ) = A 2 is 
not zero, like in the local sum-rule approach |3|, but sizeable. The (vacuum) non-locality 
parameter is the only one involved in the NLC-SR method, and it has been estimated within 
the QCD sum-rule approach from the mixed quark-gluon condensate of dimension 5: 

2 _ (g(0)V 2 g(0)) incW (gO) {ig^Gnqm _ f 0-4 ±0.1 GeV 2 [12] 

q ~ (9(0)g(0)> Hmit 2(g(0)g(0)> \ 0.5 ± 0.05 GeV 2 [13] ' ^> 

Lacking an exact knowledge of NLC of higher dimensionality, the non-locality can only 
be taken into account in the form of an ansatz. 

A simple model is provided by a Gaussian-like behavior of NLC 0, E3|, whereas A~ x 
reveals itself as the typical quark-gluon correlation length in the QCD vacuum. Hence, 
we set in coordinate space, M$(z 2 ) = (q(0)E(0, z)q(z)) ~ (qq) exp (— A 2 |2 2 |/8) ,* while the 
coordinate behavior of other NLCs looks more complicated. 



* Here, as usual, E(Q,z) = Pexp(z J* di M ^4°(t)r a ) is the Schwinger phase factor required for gauge invari- 
ance. 
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Estimates of A 2 from instanton approaches 0, 3] are somewhat larger: A 2 ~ 2/p 2 > 
0.6 GeV 2 (where p c is the characteristic size of the instanton fluctuation in the QCD vacuum: 
0.33 < p c < 0.6 Fm), while lattice calculations [IB, US 13 confirm qualitatively the Gaussian 
law for the M s (z 2 ) decay, yiel ding values close to those given in Eq. (jl.2|) . More calculational 
details are relegated to 0, [H, H3 . 

The main results of the present analysis can be summarized as follows. The moments of 
the pion DA have been re-calculated more accurately and, more importantly, their confidence 
levels have been determined. In this way, we are able to construct a whole "bunch" (a 
spectrum) of pion DAs, allowed by the NLC-SR for different values of X 2 , as noticed above. 
The important thing to be emphasized is that the range of this spectrum in the 02, 04 plane 
strongly overlaps with the 95% and also with the 68% region of Schmedding and Yakovlev 
Perhaps somewhat surprisingly, the shape of the optimum pion DA, we have determined, 
is not "dromedary" -like, as the asymptotic solution, but has two maxima (i.e., it is "camel"- 
like). Nevertheless, we stress, it is not Chernyak-Zhitnitsky (CZ) like either because in the 
endpoint region it is strongly suppressed. 



II. SPECTRUM OF PION DAs FROM NLC-SR 

Revision of pion DA moments. The NLC-SR for the DAs of the pion and the effective 
Ai + 7r'-resonance, (p n (x) and <paA x )i respectively, which appear in the "axial" channel, were 
constructed in [H, and were analyzed in 0, . Here we display only the sum rule itself: 

(U) 2 <p«{x) + (f Al ) 2 <p M (x)e" m V M2 = f ° ffrtfc s ) e - s l M2 ds + A$ s (x; M 2 ) (2.1) 

Jo 

+ A$ y (x; M 2 ) + A$ Tl (x; M 2 ) + A$ Ta (x; M 2 ) + A$ Ta (x; M 2 ) + A$ G (x; M 2 ) 

and re-estimate the moments and their error bars. The NLC contributions encoded in 
A$r(^; M 2 ) on the rhs of the SR, as well as the explanation of the individual terms are dis- 
cussed in [|&|2{J. The corrected form of the contribution for A$2\(x; M 2 ) (originating from 
the NLC (qGq)) and also the main contribution to the rhs of the SR, which is accumulated 
in the A$s(£; M 2 ) term (owing to the NLC Ms), are shown in Appendix A. 

Both the sensitivity and the stability of the NLC-SR are considerably improved relative to 
the standard SR (but nevertheless we still assume a rather conservative size of the errors of 
the order of 10%). This allows us to determine the first ten moments {£ n )-k = Jq (p 7V (x)(2x — 
l) N dx of the pion DA quite accurately. This is illustrated in Fig. QJa), where the curves for 
(£ 4 )?r are represented as functions of the Borel parameter M 2 plotted within the range of the 
stability window of the SR. The solid line (optimum moment) is determined from the NLC- 
SR (|2.ip . for so = 2.2 GeV 2 , while the broken lines give the results for thresholds varied by 
10% around this value, which turns out to be practically independent of the moment order N 
Then, the fidelity windows are: 0.5 GeV 2 < M 2 < 2.0 GeV 2 for all N = 0, 2, ... , 10. 
Notice that the range of stability and the Borel fidelity windows almost coincide with each 
other, starting for all N at M 2 w 0.6 GeV 2 . 
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FIG. 1: (a) Original SR for the moment (£ 4 )tt and (b) "daughter SR" for the inverse moment (cc _1 )^ R , both 
as functions of the Borel parameter M 2 , as obtained from the NLC-SRs. Both kinds of SRs were processed 
by including the Ai-meson with so = 2.2 GeV 2 . For the depicted range of M 2 , the fidelity and stability 
windows almost coincide. Solid lines correspond to the optimal threshold So, whereas broken lines denote 
the results for the upper (long-dashed line) and lower (short-dashed line) moment limit with a 10%- variation 
of s . 
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0.086 
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NLC SR 19} : n 
NLC SR [20] : vr 


0.131 


0.25 


0.12 


0.07 








0.131(2) 


0.25(1) 


0.110(7) 


0.054(3) 


0.031(2) 


0.0217(7) 


2.75(5) SR 


NLC this work : tt 
NLC this work : A\ + tt' 
CZ SR 22] : tt 


0.131(8) 


0.265(20) 


0.115(12) 


0.061(8) 


0.037(5) 


0.024(4) 


3.35(32) SR 


0.210(17) 


0.21(2) 


0.116(12) 


0.078(8) 


0.055(6) 


0.042(5) 


3.6(4) SR 


0.131 


0.40 


0.24 











Table 1. The moments (S, N ) m(^ 2 ) determined at ji 2 ~ 1 GeV 2 with associated errors put in paren- 
theses. Recall that CZ give all moments normalized at the lower scale /i 2 , = 0.5 GeV 2 . 

The important convolution (x~ 1 ) 7T = / — dx, appears in the perturbative calcula- 

Jo x 

tion of the 7*7 — > n° process. We construct a "daughter SR" directly for this quantity 
from Eq. (J2.1)) by integrating its rhs with the weight 1/x. Due to the smooth behavior 
of the NLC at the end points x = 0, 1, this integral is well defined, supplying us with an 
independent SR, with a rather good stability behavior of (a; _1 )^ R (M 2 ), as one sees from 
Fig. E^b). To understand this behavior, let us recall an important feature of the NLC, 
encoded in A$r(£; M 2 ). This contribution is suppressed in the vicinity of the end-points 
x = 0, 1, the range of suppression being controlled by the value of the non-locality param- 
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eter A 2 . The larger this parameter, at fixed resolution scale M 2 > A 2 , the stronger the 
suppression of the NLC contribution. Similarly, an excess of the value of (x~ 1 ) 7T over 3 
(asymptotic DA) is also controlled by the value of A 2 , becoming smaller with increasing 
A 2 .. Note, that instanton-based models demonstrate the same behavior of the DA with a 
non- locality parameter proportional to p~ 2 jl^]. 

From Table 1 one infers that the mean values of the re-calculated moments are close to 
the old ones, whereas the value of (x _1 )| R changes significantly. The reason is the corrected 
term (x; M 2 ) (in the NLC-SR), which allows us now to reconcile this quantity with the 
value obtained by calculating it with the model DAs (see below) . The new estimates of the 
error bars in Table 1 have been formed from different sources. Note that the most significant 
error (about 50%) stems from the uncertainties of the parameters of the A\ + 7r'-resonance. 

Spectrum of admissible pion DAs. Models for the pion DA, in correspondence to 
the moments in Table 1, can be constructed in different ways 0, However, it appears 
that two-parameter models, the parameters being the Gegenbauer coefficients a2 and (as 
also used in 0), enable one to fit all the moment constraints for (C, N ) n given in the Table, 
as well as to reproduce the value of (x -1 )^ within the quoted error ranged The optimum 
model DA is 

1 + a° ptl • Cl /2 (2x - 1) + af u • cf 2 (2x - 1)1 , (2.2) 



<^ pt (*) = ^ s (x) 



a° ptl = +0.188 , a4 ptl = -0.130 , with X 2 ~ 10~ 3 , (2.3) 

yielding (ar" 1 )^* = 3.174. But one can construct a whole admissible set (a spectrum), 
{fi(x; 02, cla)}, of such models by demanding that the associated moments lie inside the error 
bars, presented in Table 1. This set of DAs approximately corresponds to those parameters 
(02,04), which satisfy x 2 < 1- I* 1 this way, we obtain a "bunch" of admissible DA profiles, 
shown in Fig. El^a) in terms of dashed lines, in addition to the optimum one (thick solid 
line). For these DAs, the corresponding values of {x~ l ) Vl vary in the interval 

3.08 < (ar 1 )^ < 3.24 vs 3.03 < {x~ x ) s ^ < 3.67 . (2.4) 

Let us mention at this point that in the standard approach the value of (x _1 )(yU 2 ) cannot 
be determined from the sum rule itself due to the presence of singular 5(x)-terms. On 
the other hand, for the CZ model DA, one finds (a?~ 1 )czG u2 ) ~ 5.1 at normalization scale 
fi 2 « (0.7) 2 GeV 2 or « 4.44 at the scale fi 2 « 1 GeV 2 . Meanwhile, our estimates of the 
inverse moment in ()2.4|) . as well as the value of 02 in ()2.3j) are in a good agreement with the 
estimates (ar 1 )^ GeV 2 ) = 3.3 ± 0.3, and a 2 (l GeV 2 ) = 0.1 ± 0.1 obtained in from an 
analysis of electromagnetic pion form factor. 

The above results have been obtained within the framework of NLCs, employing a single 
non-locality parameter, viz., A 2 , fixed at the so-called "standard" value of 0.4 GeV 2 [12]. 
To examine the role of this parameter in determining the structure of the model DAs, we 
constructed another set of amplitudes, <f2(x), fixing it at the higher, yet admissible, value 
A 2 = 0.5 GeV 2 13]. Again, an "optimum" DA, Lf2 Pt (x), was determined with coefficients 
o P t2 _ _|_q a^ 2 = —0.091 with a profile similar to that of <£>° pt (x), as Fig. |^b) reveals. 



a 



t Note that SY assume that a„>4 are small, while in our approach these coefficients have been calculated 
up to N = 10 and found to be numerically negligible. 
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FIG. 2: Graphical representation of <pi(x) (part (a)) and <fi2(x) (part (b)). The thick solid lines in both 
plots denote <£>° p 2 (x), i.e., the best fit to the determined values of the moments (see Table 1), whereas dashed 
lines illustrate admissible options with \ 2 < 1- 

Note, however, that now the envelope of the "bunch" becomes quite close to the asymp- 
totic DA, i.e., "dromedary-like (convex) at x = 1/2 and correspondingly it exhibits less 
suppression in the endpoint region x = 0,1. We note in passing that a suppression of the 
endpoint region as strong as possible (for a discussion we refer to |25j ) is important in order 
to improve the self-consistency of perturbative QCD in convoluting the pion DA with the 
specific hard-scattering amplitude of particular exclusive processes. As before, the values 
of the inverse moment (x -1 )^, extracted from the new SR, and those following from the 
model f2(x), are mutually consistent: 

3.05 < (x~% 2 < 3.22 vs 2.87 < (x' 1 )^ < 3.51 . (2.5) 

It is worth remarking here that this "camel-like" structure of the "best fit" (alias, the 
optimum) DAs for both considered values of the non- locality parameter A 2 ,, displayed in 
Fig. |2ta,b) can actually be traced back to a definite origin. Crudely speaking, this shape 
structure is the net result of the interplay between the perturbative contribution and the 
non-perturbative term A&si x ) that dominates the rhs of the SR in Eq. (|2.1|) . The fact that 
the function A<&s{ x ) is n °t singular in x and has a dip at the central point of the interval 
[0, 1] is also reflected in the shapes of the DAs. 



III. QCD SR RESULTS VS CLEO DATA 

Overlap of regions in plane. Let us now draw in Fig. a) the regions of admissible 
pairs (a 2 , a 4 ), determined in correspondence with the error bars of the moments, in order 
to compare them with the experimental constraints supplied by the recent high-precision 
CLEO data [4] . The latter are processed in theplot in Fig. E^b) in the form of confidence 
regions, extracted by Schmedding and Yakovlev [5J] using a NLO light-cone QCD SR analysis 
that (approximately) corresponds to an average normalization point of /i = 2.4 GeV. 

The region, corresponding to the bunch of DAs, determined at the "standard" value 
A 2 = 0.4 GeV 2 , and evolved to /i = 2.4 GeV in Fig.|2fa), is enclosed in Fig. |Sfa) by a slanted 
rectangular area. Superimposing this area with the "ellipses" pertaining to a 68% (solid 
line) and a 95% (dashed line) confidence region, [5] , we show in Fig. |3fb) to which extent 
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FIG. 3: (a) The parameter space of (02, cla) pairs (enclosed by a thick solid line), corresponding to the 
allowed values of the second and fourth Gegenbauer coefficients, calculated within the NLC-SR approach for 
three different values of A 2 at [i = 1 GeV. The optimum DA is denoted by a black dot. (b) shows how our 
estimate of the confidence region with A 2 — 0.4 GeV 2 at fi = 2.4 GeV overlaps with those displayed in Fig. 6 
of Ref. 5]. Bold dots in plot (b) mark the parameter pairs for the asymptotic (big black dot), Chcrnyak- 
Zhitnitsky (full square) and our optimum (with respect to the SR) DA (black dot). Contour lines show 
68% (solid line) and 95% (dashed line) confidential regions, extracted by Schmedding and Yakovlev from 
3000 randomly chosen sets of CLEO data 4jj . The small axis of the "ellipses" is generated by experimental- 
statistical uncertainties, whereas the large one results from theoretical-systematical uncertainties. 



our estimates for the pion DAs agree with the CLEO data. In particular, the central point 
(a^, fl4 Pt ) of the slanted rectangle lies just on the 68%-region and the lower corner of it 
"touches" the central point ( ) of the SY plot: 

a* = 0.190, a* = -0.14 , (3.1) 
a°, ptl = 0.139, a^ 1 = -0.082 . (3.2) 

Given the errors determined by SY, 

a 2 = a* 2 ± 0.04(stat) ± 0.09(syst), a 4 = a* A ± 0.03(stat) =f 0.09(syst) , (3.3) 

one realizes that our theoretical estimates are quite compatible with their experimental 
constraints. The region enclosed by the slanted rectangle bounded by the long-dashed line 
along the "diagonal" (see below) in Fig. E(a), corresponds to the bunch of DAs displayed 
in Fig. Efb). Though still within the 95% confidence region of SY, it is, however, mostly 
outside the central 68% region. Finally, the third slanted rectangle limited by the short- 
dashed contour, and shifted along the "diagonal" to the upper left corner of the figure in 
Fig. Of a), corresponds to a trial NLC-DA with A 2 = 0.6 GeV 2 . This value falls actually 
outside the standard QCD NLC-SR bounds in Eq. (|1.2J) for A 2 . Remarkably, the image of 
this region in Fig. 0(b) (not displayed) would lie completely outside the central region as a 
whole. Therefore, we conclude that the CLEO data does not prefer the value A 2 = 0.6 GeV 2 , 
in full agreement with previous QCD SR estimates."'" 



* A very interesting detail in Fig. E(a) is that the intersection set of all three admissible regions represented 
by the slanted rectangles is located around (and contains) the optimum point for the value A 2 = 0.5 GeV 2 , 
though its image in Fig. |3f b) lies outside, but close to the upper boundary of the 68% SY region. 
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Role of the "diagonal". By measuring the 7 + 7* — > it form factor, Fjj* w , one mainly 
probes the pion inverse moment {x~\) n = 3(1 + a 2 + a 4 + . . .). To illustrate this, let us apply 



the leading twist NLA expression [26] to the Q 2 F, 



3 

—Q 2 (F 77 .^ « C ® <p v ) = I = 3 [1 + a 2 + a 4 + . . . + a s (A + a 2 A 2 + a 4 A 4 + ...)] . (3.4) 

Here C is a coefficient function in the factorized amplitude for the process at NLO; A n 
gives the size of the leading radiative corrections to the contribution of the Nth Gegenbauer 
eigenfunction entering the expansion of ^p n (x). It should be noted that this standard fac- 
torization formula is not well established when one of the photons in the process is on-shell. 
Nevertheless, for demonstration purposes this approach is good enough. The correspondence 

/(/i 2 ) = (x-%{v 2 ) + 3« s (/i 2 )A (3.5) 

is slightly smeared by radiative corrections due to higher eigenfunctions in Eq. (J3.4)) . There- 
fore, constraints imposed by F^*^ measurements are in reality constraints for the combina- 
tion a 2 + a 4 + . . ., modulo contributions from radiative corrections and higher twists. Due to 
this reason, the admissible region for the desired points (02,04) in 0, (Fig. 0(b)), is, so to 
speak, stretched along a "diagonal" defined by a 2 + 0,4 = const. The value of this constant 
can be easily extracted from the data with a reasonable accuracy,^ using Eqs. (|3.2[) . (j3.3J) . 
to obtain 

a 2 + a 4 = 0.05 ±0.07 (3.6) 
(large systematic errors cancel!) that should be compared with our value 

a 2 + a 4 = 0.056 ± 0.033 . (3.7) 

Hence, we see that an "experimental" quantity (Eq. ([3.6)1 ) expressively agrees with those 
obtained in the NLC-SR (cf. Eq. (|3.7|) ). On the other hand, the direct estimate for (x^ 1 )^ 
from this SR (at /j 2 ~ 1 GeV 2 in Table 1) is also in agreement with both of these estimates, 
(£~ 1 )tt/3 — 1 ~ 0.1 ±0.1. The evolution of this quantity to the scale n 2 = (2.4 GeV) 2 should 
eventually make the agreement more pronounced. To complete these considerations, let us 
mention that the empirical SY estimate, introduced in [Hf for a "deformed diagonal" in the 
(a 2 , c^-plane 

a 2 + 0.6a 4 = 0.11 ±0.03 , (3.8) 
is also satisfied by the values of the central point of our admissible region, namely, 

a 2 + 0.6a 4 = 0.09 ± 0.039 . (3.9) 

We emphasize that the radiative corrections in Eqs. (|3.4j) . (|3.5J) are particularly important 
because it turns out that they dominate over the "non-perturbative" a 2 , a 4 contribution. To 
make this point more apparent, let us consider the asymptotic DA in NLO approximation 
[I?} (see Eqs. ()B.2|) . ()b!3|) in Appendix B). Then, we find that this eigenfunction supplies 
the main contribution to / and a s A : 



An 



-5C F + [ 2C F - ^ 



2^ 



-7 -0.17. (3.10) 

47T 



At the same time, the theoretical-systematical errors to every pair of values of 02,04 appear to be very 
large, see Eq. (|3.3|l . 
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The first contribution, — 5CV, originates from the coefficient function (see |28|). while the 
second one, (2CV — &o/3), from the modification of the asymptotic DA at NLO j22|. The 
final estimate, —0.17, corresponds to the value a s /(47r) = 0.024 at the normalization 
point n = 2.4 GeV. We see that the estimate for a s A is indeed 3 times larger than those 
for a 2 + a 4 , considered above, and for that reason it controls the level of accuracy in the 
determination of the latter.^ The twist-4 contribution to the form factor is also important 
(see (il), being of the same order as a 2 + a 4 but having the opposite sign. The overall 
radiative correction to the contribution of still higher eigenfunctions is also negative and 
small in comparison to the estimate in Eq. (j3.10j) . 

^Jil ■ (a 2 A 2 + a 4 A 4 ) « 0.024 • (-0.58) = -0.014 . 
Atx 

To further improve the theoretical accuracy of extracting a 2 and a 4 (from experimental data), 
one should first obtain an estimate for C (g> y? as in next to NLO approximation because this 
contribution will certainly dominate in that order. 



IV. COMPARISON WITH RESULTS FROM DIFFERENT SOURCES. 

In this section we compare the results of our analysis with new evidences for the form 
of the pion DA, following from hadron phenomenology, lattice simulations and instanton- 
based models. However, keep in mind that the information provided from these sources 
bears uncertainties the size of which is difficult to estimate. 

DA from 7r into di-jets via diffractive dissociation. It is tempting to compare the 
results obtained and discussed above with the intriguing interpretation of experimental data 
on di-jet production in pion-nucleon interactions |6j. This process was suggested already 
in (3^ as a tool to measure the shape of the pion DA directly from the data via diffractive 
dissociation into two jets. Such an experiment was carried out very recently by the Fermilab 
E791 Collaboration [6( by measuring the transverse momentum distribution of the diffractive 
di-jets in two windows of k t : 1.25 < k t < 1.5 GeV and 1.5 < k t < 2.5 GeV . We have found 
that their experimental deduction for the shape of the pion DA in the lower window can be 
well approximated by a model with two Gegenbauer coefficients to read 



<f%(x, Q 2 ) = ^(x) [1 + 4\Q 2 ) ■ Cl /2 (2x - 1) + af(Q 2 ) ■ C 3 A /2 (2x - l)j , (4.1) 

with x 2 ~ 10~ 2 and af = 0.121, af = 0.012 at Q 2 « 8 GeV 2 . Despite the large "diagonal" 
value af(fi 2 ) + af(/j, 2 ) « 0.143 (at fi = 2.4 GeV), resulting from this model, it is still 



" An analysis in leading approximation, performed in 30] for the phenomenological value I sa 2.4, leads to 
the conclusion that / w (x^ 1 )^, so that the associated DA is forced to be narrower than the asymptotic 
one. But in NLO this estimate transforms into Eqs. (|3.4|) . (|3.5f) and therefore one can conclude that 
the value of (x -1 )^ is larger than 3, with the pion DA being broader than the asymptotic one. Indeed, 
Fig. EHb) demonstrates exactly this fact in detail. In the same context we mention that a narrower DA 
than the asymptotic one was also used in the calculation of i 7 ! n - 7 in yielding a prediction somewhat 
below the CLEO data and leading to the suggestion that the shape of the "true" pion DA must be broader 
than that of the asymptotic one. 
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consistent with the SY constraints, as one realizes by comparing it with the values given 
in Eqs. (13. 6J) . ()3.7|) . The symbol "star", associated with ^(/i 2 ), is located in Fig. Hfb) 
on the upper boundary of the SY 95% region and above the diagonal. Fig. HJa) shows 
the profile of the E791-fft (thick solid line) in comparison with our "optimum" DAs for 
A 2 = 0.4 (dashed line) and A 2 = 0.5 GeV 2 (solid line), both evolved to Q 2 « 8 GeV 2 . 
Taking into account that the accuracy of the data is not better then 10%, one is tempted to 
conclude that the DA shape from the di-jet data is in fact not far away from the yvbunch. 
Note that the authors of the second paper in [32| have discussed this treatment of di-jet 
production data and concluded that the interpretation of them with respect to the pion DA 
is yet questionable. Indeed, inspection of their Fig. 3 p reveals that the results (data) for 
x < 0.1 and x > 0.85 — 0.9 are different, exhibiting the limitations of the accuracy of their 
fit in the endpoint regions. 

Transverse lattice. Quite recently, Dalley |7[ has used the data points in a transverse- 
lattice calculation to deduce the shape of the pion DA and approximated it by the following 
analytic form 



V lat (x, /i 2 ) = </? as (x) 1 + 0.133 • Cl /2 (2x - 1) (4.2) 

at /iQ ~ 1 GeV 2 . The shape parameters 02,04 of this DA are inside the admissible 95% 
region in Fig. EJ^b) (black triangle) and somewhat closer to the diagonal than the fit of 
E791. Because of several approximations in this calculation, the quality of the deduced fit 
is unknown and as a result the errors in the coefficients are presumably large. Hence, we 
may regard this model DA as being only qualitatively consistent with the y^-bunch. 

Two things are worth to be remarked. First, this new model fit is in conflict with the 
DA shape predicted in |33[ and more recently in [34| using the same Hamiltonian and 
Fock space truncation, but employing continuous basis functions rather than the discrete 
light-cone quantization applied by Dalley. Second, if the shape obtained is very close to 
the asymptotic one, then it is again ruled out by the CLEO data, as it is outside the SY 
confidence regions. For the sake of completeness, we mention in this context that other 
lattice attempts 3^, 3(| yield unrealistically small values of the moment (t; 2 )-*- 



Instanton-induced models. Instanton-induced models, worked out in [37, 38], usually 
lead to DA shapes close to the asymptotic one. As a consequence, such DAs are definitely 
ruled out by the SY constraints shown in Fig. EJb). However, these DAs satisfy "a softer 
constraint" , also derived by SY in (see Fig. 7 there). In fact, the confidence regions derived 
this way, are somewhat enlarged with the effect that the asymptotic DA is just "inside" the 
larger region and positioned at the upper end of the diagonal. Note that instanton-induced 
models prefer a higher value of the vacuum quark virtuality: A 2 > 0.6 GeV 2 . For such a value 
of the parameter A 2 (unrealistically large for QCD SRs), the condensate contribution on the 
rhs of Eq. (|2.1j) becomes small in comparison to the perturbative one. As a consequence of 
this, the corresponding DA appears to be closer to the asymptotic one, as we have already 
indicated (cf. Fig. Eta), where one sees that the center of the slanted rectangle for this A 2 
value is much closer to the o 4 axis and hence to the asymptotic DA.) 

A separate case of an instanton-induced pion DA, consistent with the SY constraints, has 
been derived by Praszalowicz [39]. We included this model in Fig. Efb) by the label "full 
diamond" . 
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FIG. 4: (a) Selected profiles of ip w (x,/i 2 — 8 GeV 2 ). The short (long) dashed line corresponds to our 
"optimum" DA model, ^° pt (x) at A 2 = 0.4 GeV 2 (<p?\x) at 0.5 GeV 2 ) and the solid line to the fit to the 
model deduced from the experimental data on di-jet production in pion-nucleon interactions [g]. (b) This 
plot is obtained by inserting into Fig. 6 of Ref.0 the confidence region determined for the (^2-bunch (at 
A 2 = 0.5 GeV 2 ) and confined in the slanted rectangle. The E791 Collaboration fitQ, Eq. I|4.1|) . is denoted 
by a black star, whereas the Dalley model, Eq. I|4.2|l . is marked by a black triangle. All model DAs shown 
are evolved to \i = 2.4 GeV. 



V. CONCLUSIONS 

1. We performed an accurate analysis of QCD SR with NLCs for the pion DA and 
obtained admissible sets of them with respect to the SR constraints for different values of the 
non-locality parameter A 2 ,. Confronting these determined DA spectra with the constraints 
obtained by SY in [^] from the high-precision CLEO experimental data [4], we showed that 
they satisfy them quantitatively (cf. Fig. Eland Fig. El^b)). 

2. We found that the predictions for the profile of the pion DA deduced from the di-jet 
production (E791 Collaboration) and those extracted from the simulation on a transverse 
lattice 0] are consistent only with the 95% SY constraints (2<j-deviation criterium), but not 
with the 68% SY constraints (la- deviation criterium). Moreover, they are located near 
the boundary of the 2cr-region obtained with the standard error estimation. As shown in 
Fig-IHb), the di-jet prediction is not too far from the yvbunch region (with A 2 = 0.5 GeV 2 ), 
though a quantitative comparison calls for a higher accuracy of the data processing. On the 
other hand, the lattice result qualitatively also agrees with the set of our DAs (the y? 2 -bunch) 
(see Fig. Efb)), but a more detailed comparison is obscured by the unknown quality of this 
method. 

3. The existing predictions from instanton-induced models are too close to the asymptotic 
DA and therefore they are only consistent with the "softened" version of the SY constraints 
(see for details in 0), except for the new model by Praszalowicz |39[, which is just outside 
the confidence region of the yvbunch and on the boundary of the 95%-region of SY. 
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APPENDIX A: EXPRESSIONS FOR NONLOCAL CONTRIBUTIONS TO SR 



The form of contributions of NLCs to the OPE on the rhs of Eq. (j2.1|) depends on the 
modelling of NLCs. At the same time, the final results of the evaluation of the SRs show 
stability against variations of this modelling, provided the scale of the average vacuum quark 



virtuality A 2 is fixed. Here, we used the model (delta-ansatz) suggested in (3, llCj and used 



21 



This model leads to a Gaussian decay of the scalar quark condensate 
Ms(z 2 ). For the four-quark condensate, the factorization ansatz is applied to reduce its 
contribution to a pair of scalar condensates. In the NLC approach this may lead to an 
overestimate of the four-quark condensate contribution, A$s (x\ M 2 ), because it evidently 

A 2 /(2M 2 ), A = 1 - A: 



extensively in |2 

v2 



neglects the correlation between these pairs. Below, A 



A$ s (x; M 2 ) 



A s 18 
M 4 AA 2 



Instead of (A. 5) in |20j], we have 



9 (x > A > x) x [x + (A — x) In (x)} + (x — > x) + 

-6(1 > A)6 (A > x > A) [A + (A - 2xx) ln(A)] } . (A.l) 



A$ Tl (x; M 2 ) 



■^0(1 > 2A) i[8{x - 2 A) - S(x - A)] (i - 2 



e(x > 

A 



x-2A 
AA 



+ (x 



X) 



9{2A > x)- 

(A.2) 



Here, A s = (87r/81) (\/aIq(0)q(0)) 2 , whereas for the quark and gluon condensates we use 
the standard estimates (v^(O)g(O)) w (-0.238 GeV) 3 , (a s GG) /12n w 0.001 GeV 4 



APPENDIX B: RADIATIVE CORRECTIONS 

The radiative corrections to the correlator and to the pion DA are 



p^\x,s) = 3xx< 


I- 


a s 


c F \ 


^ s LA (x) = 6xx < 


!■♦ 


a s 




(x- 1 )^ = 31- 


f a s 





(B.2) 
(B.3) 



Here, bo is the first /3-function coefficient, bo = (11/3)Ca — (4/3)TrA^/; a s = a s (/i 2 )/47r, 
A^° w 0.465 GeV, af LO (m 2 ) w 0.358, see, e.g., Q. We perform a 2-loop evolution of the 
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pion DA following the approach presented in the first paper in [26| and using the so-called 
"optimized MS-scheme" 41 . 
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Erratum 



We used a wrong formula, Eq. (B.l), for the radiative corrections to the spectral density 
of the axial-axial correlator. The correct one should read as follows (x = 1 — x): 



p pert (x,s) = 3 xx 



7Y 



1 + a s C F ( 5 - '— + ln 2 (x/a;) 



1 



2tt 2 



(2.4) 



We reanalyzed our sum rules using the correct expression and obtained new parameters for 
the phenomenological spectral density in the axial channel: 



(2.5) 



where Ai in fact represents a mixture of the real A\ and n' mesons. The optimum stability 
of the analyzed sum rules is achieved for s = 2.25 GeV 2 and m 2 Ai = 1.616 GeV 2 . The 
error induces mainly a shift in the determination of the decay constants of the pion and its 
excitations; it does practically not influence the values of the extracted moments of the pion 
distribution amplitude, as one can see from Table 1, where we collect the values of these 
quantities in comparison with the published ones. 



U (GeV) 


N = 2 


N = 4 


N = 6 


N = 8 


N = 10 


(O 


0.131(8) 


0.265(20) 


0.115(12) 


0.061(8) 


0.037(5) 


0.024(4) 


3.35(32) SR 


0.137(8) 


0.266(20) 


0.115(11) 


0.060(7) 


0.036(5) 


0.025(4) 


3.35(30) SR 



Table 1. The moments (^ n )tt{p- 2 ) determined at fi 2 ~ 1.35 GeV 2 with associated errors put in 
parentheses. The old values are given in the first row and the new ones in the second row. 



Table 2 represents the changes for the effective Ai-meson decay constants and moments: 



f Al (GeV) 


N = 2 


N = 4 


N = 6 


N = 8 


N = 10 


(x- 1 ) 


0.210(17) 


0.21(2) 


0.116(12) 


0.078(8) 


0.055(6) 


0.042(5) 


3.6(4) SR 


0.221(20) 


0.21(2) 


0.113(12) 


0.076(8) 


0.055(6) 


0.040(5) 


3.5(4) SR 



Table 2. The moments (£ N } Al (fi 2 ) determined at ^ 2 - 1.35 GeV 2 . 



The old and new parameter regions of (02, 04) pairs, corresponding to the allowed values 
of the second and fourth Gegenbauer coefficients, calculated with the corrected sum rules 
in comparison with the published ones for the value A 2 , = 0.4 GeV 2 at /i 2 ~ 1.35 GeV 2 are 
given in Fig. 1. 

Given that the two (a 2 , a A ) areas almost coincide and the BMS distribution amplitude is 
only infinitesimally shifted none of the results of our analysis is affected. 
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FIG. 5: Corrected (02, 0,4) area (blue solid line) in comparison with the published one (red dotted line). 
The blue cross marks the new position of the BMS distribution amplitude, while the red circle denotes the 
old one. 



16 



